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T^, ' Abstract. Bethe-Salpeter and light-front bound state equations for three 

scalar particles interacting by scalar exchange-bosons are solved in ladder 
^ , truncation. In contrast to two-body systems, the three-body binding energies 

' obtained in these two approaches differ significantly from each other: the 

■ ladder kernel in light-front dynamics underbinds by approximately a factor of 

two compared to the ladder Bethe-Salpeter equation. By taking into account 
three-body forces in the light-front approach, generated by two exchange- 
bosons in flight, we find that most of this difference disappears; for small 

OO , exchange masses, the obtained binding energies coincide with each other. 

O 

^ . 1 Introduction 

5t , Relativistic two-body systems are often treated using the ladder truncation in the 

Bethe-Salpeter (BS) framework [1] or in light-front dynamics (LFD) [2l|3]. These 
two approaches give results very close to each other for the binding energy [1] 
and also for electromagnetic form factors [5] in the spinless case. Binding energies 
of a two-fermion system were calculated both in the BS approach in [6] and in 
LFD [7] and again the results turned out to be rather close to each other. 

In these cases, the relativistic calculations using either the BS approach or 
LFD lead to binding energies that are quite different from the non-relativistic 
binding energy, obtained by solving the Schrodinger equation. The relativistic 
effects are more important for larger binding energy. In addition, they are also 
important for small binding energy, if the mass of the exchange particle is large 
enough, that is, of the order of the constituent mass [4]. 

Here we extend the same methods, the BS approach and LFD, to scalar three- 
body bound states, using a one-boson exchange kernel. Eventually, the goal is to 
extend these methods to bound states of fermions: e.g. using a meson-exchange 
model in the case of 3-nucleon bound states (triton, ^He), and baryons as bound 
states of three (non-perturbatively dressed) quarks interacting via gluons. For 
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(a) 




(b) 

Figure 1. (a) Top: Feynman ladder graph with two exchanges, (b) Bottom: One of six time- 
ordered graphs, generated by the ladder Feynman graph (a). 



simplicity however, and as a first step, we restrict ourselves here to spinless 
systems. Prehminary results for the three-body BS equation with a one-boson 
exchange kernel have been presented in [8]. As far as we know, the LFD three- 
body bound state equations with one-boson exchange kernel has never been 
solved for such systems. Previously, the LFD equation was solved in [9l [10] for 
zero-range interaction. Its solution was also found in the relativistic quantum- 
mechanical approach for scattering states [11], with a phenomenological mass 
operator. 

It is important to note that the actual ladder kernels in the BS approach 
and in LFD are not identical, nor are they given by the same graphs. In the BS 
ladder kernel, there is no notion of time-ordering in the diagrams. On the other 
hand, the kernel for the LFD equation is given by the time-ordered graphs in 
the light-front (LF) time. 

Thus, the second iteration in the two-body ladder BS equation (the second- 
order Feynman box graph. Fig. [H^a)), when represented as a set of time-ordered 
graphs, turns into six LF time-ordered graphs, including two so-called "stretched 
boxes" with two exchange particles in the intermediate state. One of such 
stretched boxes is shown in Fig. Wijo)- These stretched boxes with two (and 
more) exchange particles in the intermediate state are implicitly included in the 
BS equation, but they are not generated by iteration of the LF ladder kernel; 
therefore, they are omitted in the corresponding LFD bound state equation. In 
principle, this will cause a difference between the BS and LFD results. However, 
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direct examination of the stretched boxes [12] shows that they are small. Being 
explicitly included in the LF kernel, they result in a small correction to the bind- 
ing energy [13]. This explains, why the BS and LFD two-body binding energies 
are very close to each other, despite the fact that the kernels are not identical. 

In the three-body problem, the situation is quite different: even with a sim- 
ple one-boson exchange kernel, new interesting and important features emerge, 
highlighting differences between the BS approach and LFD. To be specific, in a 
three-body problem the difference between LFD and time-ordered iterated BS 
kernels cannot be reduced to stretched boxes: new LFD diagrams appear, which 
are shown in Fig. [2j Like stretched boxes, these diagrams cannot be constructed 
by iterations of the two-body LF kernel. In contrast to the stretched boxes, they 
do not contain any loops. These graphs generate three-body forces which are not 
taken into account in the three-body LFD equation with ladder kernel. On the 
other hand, they are implicitly included in the three-body BS equation. 

Three-body nucleon forces have been discussed for several decades. In recent 
years, there has been an increased interest in 3-nucleon forces due to high preci- 
sion few-body calculations. The need for three-body forces results from nuclear 
phenomenology (see for review Ref. [H]). In particular, three-body forces remove 
the underbinding of light nuclei found with local two-body forces. With nonlocal 
two-body potentials there seems to be less of a need for three-body forces from a 
phenomenological point of view [15]. The diagrams depicted in Fig. [2] represent 
a part of three-body forces for nucleons. Another part arises due to excitation 
of intermediate isobars (most notably Z\-excitations). The first-order relativistic 
corrections due to the first eight graphs of Fig. [2l was first found in Ref. [16] . Its 
contribution to the triton binding energy was calculated in Ref. [17] . In Ref. [18] 
it was shown that this relativistic correction cancels with the corresponding rel- 
ativistic correction to the second iteration of the one-boson exchange. Due to 
this cancellation, the sum of these two corrections does not contribute in the 
Schrodinger equation. However, in a truly relativistic framework, the full graphs 
Fig. [2] (not a first-order relativistic correction) should be taken into account. 

Both the BS approach and LFD provide us with a relativistic framework to 
study relativistic three-body bound states. With modern computer resources, we 
can now solve these bound state equations in ladder truncation without further 
approximations (at least for spinless systems). These calculations serve as a the- 
oretical laboratory where we can study the importance of three-body forces of 
relativistic origin, like the ones generated by the diagrams in Fig. [2j By contrast- 
ing results from the BS approach and LFD we can elucidate their contribution 
(if any) to the binding energy. Note that in the scalar model under consideration 
here, we have no analog of the three-body forces originating from isobar exci- 
tations, so we cannot asses the relative importance of three-body forces due to 
isobar excitations; nor do we have any means to asses the role of intrinsic three- 
body forces such as NNN — > NNN point interactions, or, in QCD, three-body 
forces due to the triple gluon vertex. 

The aim of this paper is two-fold: (i) First of all, we will solve, for the first 
time, the three-body bound state BS and LFD equations with a one-boson ex- 
change kernel and compare the corresponding binding energies. We shall see 
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that, in contrast to the two-body case, the binding energies found in these two 
approaches do not coincide with each. The LFD binding energy is significantly 
smaller than the BS binding energy (i.e. the LF ladder kernel underbinds com- 
pared to the BS ladder kernel), (ii) Then, in order to identify explicitly the 
origin of this difference, we calculate perturbatively the contribution to the LFD 
binding energy from the three-body forces depicted in Fig. [2j It turns out that 
this correction eliminates most of the difference. In this way, we see a clear and 
undoubted manifestation of the three-body forces and establish that they are 
responsible for the underbinding of the LF ladder kernel. 

The plan of the paper is the following: In Sec. [2] we present the three-body 
BS equation, and we give a brief outline of our numerical procedure. The corre- 
sponding equation in LFD is derived in SecEJ the three-body kernel determined 
by the graphs Fig. [2] and its perturbative contribution to the binding energy are 
also given in this section. Next, we present our numerical calculations in Sec. HI 
where we contrast the results in LFD with those obtained from the BS equa- 
tion. Finally, Sec. [5] contains concluding remarks. Additional details regarding 
the relativistic LF Jacobi variables, used in our calculations, can be found in 
Appendix A. 



2 Three-body Bethe Salpeter equation 

A three-body bound state with total four-momentum P can be described by 
solution of the three-body bound state equatioE0 



, (27r)4y (27r)4y (27r)4 ' ^ 
X K{pi,p2,P3;h,k2,k3;P) A{ki) A{k2) Aik^) r{ki,k2,k-y, P) 

at = — M^, where M is the three-body bound state mass. Here, K is the three- 
body scattering kernel, which can be decomposed into three 2-body kernels Kij 
and an intrinsic 3-body kernel i^i23, see Fig. [3l A are the (dressed) propagators 
for the constituent particles. Momentum conservation dictates Pi + P2 + Ps = 
ki + k2 + k^ = P, which is enforced (on the right-hand-side) by the (5-function. 
In the absence of intrinsic three-body kernels {K123 in Fig. [2D, the three-body 
bound state equation reduces to 

/d^k 
j^Ki2{puP2;ki,k2;P) A{ki) A{k2) r{ki,k2,Pz]P) 

+ j -^K23{p2,P3]k2,k3-,P) A{k2) Aih) r{pi,k2,k3;P) 

+ J ■^:^K3i{p3,Pi]k3,ki-,P) Aiks) Aih) rCki,P2,k3-,P), 

(2) 

where ki = pi + k and kj = pj — k and with Kij the same two-body kernels as in 
the two-body BS equation. 



""^We use Euclidean metric in this section, so all 4-dimensional integrations are entirely Euclidean, 
and the on-shell condition for the bound state implies = —M^. 
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Figure 3. Three-body bound state equation with two-body and three-body forces. 



2.1 Scalar three-body BS equation 

For simplicity, we consider a scalar field theory, with interaction —gcp'^ip, and 
a three-body bound state consisting of three particles (p, interacting via the 
exchange of a 93 boson. We take equal mass m for the constituents, and assume 
that they are all distinguishable; the mass of the exchange boson is fi. In the 
limit /i = this becomes the well-known Wick-Cutkosky model |19j . 

We use the ladder truncation for the three-body bound state equation. Fur- 
thermore, we do not take into account any self-energies (i.e. we use bare propa- 
gators). With these truncations, the three-body BS equation reduces to (see also 
Fig.E]) 

r{pi,P2,P3;P) = 9^Y.j (04 ^(^) ^(^1) ^(^2) r{h,k2,P3;P), (3) 

where ki = pi + k, kj = pj — k, and the sum runs (cyclic) over i,j = 1, 2, 3, i / j. 
The propagators are: scalar constituent propagator 



A{p) 

and scalar exchange propagator 

D{k) 



1 



p2 _|_ jy^I 



fc2 + ^2 



The momenta satisfy momentum conservation, so the three momenta pi are 
not independent. The bound state amplitude F is actually function of only two 
(relative) independent momenta, and one can write it as a function of only two 
independent 4- vectors 

r(.pi,P2,P3;P) = r{p,q;P). 



One can make different choices for these two relative momenta p and q; the 
physics should be independent of these choices. The most natural choice are the 
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Jacobi variables in momentum space 

Pi = iP + q + p, 

P2 = iP + q-p, 
P3 = lP-2q. 

However, this choice is not unique, and we have used other choices as well; within 
the estimated numerical errors, our results are independent of this choice (as it 
should be). 

In this notation, the ladder BS equation for r{p, q; P) becomes 



r{p,q;P) = g'J ^Dik) (4) 

X \A{lP + q+p + k) A{^P + q-p-k) r{p + k,q;P) 

+ A{lP + q-p + k) A{lP-2q-k) r{j)-\k,q+\k-P) 

+ A{\P-2q + k) A{\P + q + p-k) r{p-\k,q-\k]P) 

The momentum of the exchanged particle, namely the integration variable k, is 
the same in all three terms. 

2.2 Numerical implementation 

As mentioned, the bound state amplitude r" is a function of two independent four- 
vectors p and q, but the actual independent variables are: two radial variables, 
and g^, and three angles, 6p, 9q, and cppg. In the rest-frame of the bound state, 
P^ = [i M, 0, 0, 0], we use the notation 

Pti = p[cos(6'p),sin(6'p),0,0] 

qf, = g[cos(6'g),sin((9g)cos(0pg),sin(6'g)sin(^pg),O] 
and we have the 4-dimensional integration measure 

/foc fir fTT f'2n 

d'^k= k^ dk sm^{e)de I sin / d-f . 
Jo Jo Jo Jo 

All three angular integrations have to be done numerically, in addition to the 
radial integration. Furthermore, we need to do some interpolation on the function 
r(p, q; P) we are solving for, because we have expressions like r(p—^k, g± ^k; P) 
under the integral. Thus we have a four-dimensional integral equation for the 
bound state amplitude which is a function of five independent variables, given a 
fixed (external) P^ = — M^, with a nine-dimensional kernel. For a given choice 
of P^ = — (or of the binding energy E^, = 3m — M), we can solve this (four- 
dimensional) eigenvalue equation for the corresponding coupling constant g, or 
rather, for the dimensionless constant a = (/^/(IGvrm^); the (five-dimensional) 
eigenvector is the BS amplitude, r{p'^,q'^,9p,9y,(/)py; P'^). 

We solve this numerically by straightforward discretization of the internal and 
external variables; we do not make use of any expansion of P in partial waves, nor 
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dimensionless coupling a = g /(167rm ) 



Figure 4. (Color online) Binding energies obtained by solution of the BS equation in ladder 
truncation for three different values of the exchange- boson mass: /j, = 0.01, fj, — 0.5, and /i — 1.5 
(units set by m = 1). Solid curves are the three-body BS result; dashed curves are the two-body 
BS results. 

any other set of basis functions. No matter how we re-arrange the internal and 
external variables, we have to do some interpolation on the function F we want 
to solve for; that means that there is nothing to be gained by taking the internal 
and external grids to be the same. After a suitable transformation g-^ — > x, with 
< < cxD, and < x < 1, we split the radial interval into sub-intervals, and 
use 2-point Gaussian integration inside each sub-interval; for the external radial 
grid we use the A'" break-points (with the boundary condition that for x = 1, i.e. 

—> CO the BS amplitude vanishes). For the angular integrations we use suitable 
Gaussian integration measures: Gauss-Chebyshev for cos0, Gauss-Legendre for 
cos(j) and 7. For the external angular variables O<0p<7r,O<0g<7r, and 
< (t)pq < vr and we use equidistant grids in cos Op^q and cosippg. 

The ground state solution (that is all we are interested in right now) is 
symmetric under exchange of any of the three constituents. That means that 
r{p, q; P) = r{—p, q; P), i.e. the solution we are looking for is symmetric under 
the combined transformation 

COs(6'p) -COs{9p), COs((/)pg) -COs{(j)pq). 

Using this symmetry property saves a factor of two in the computational effort. 
Since we are only interested in the largest eigenvalue, and its corresponding 
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eigenvector (i.e. in the ground state), we can solve the eigenvalue equation by 
iteration (also called power method). Typically, we need of the order of 25 to 40 
iterations in order to obtain a stable solution. Our results are shown in Fig. [J] 
for three different values of the mass of the exchange boson. The mass m of 
the constituents is set to unity, so the maximal binding energy for a three- 
body system is three. For comparison, we also show the binding energy of the 
two-body system (also in ladder truncation), which clearly shows that the three- 
body system is stronger bound than the corresponding two-body system; i.e. the 
solutions we find are indeed three-body bound states. 



3 Three-body bound state equation in LFD 

The most simple way to derive the LF equation for spinless particles and calculate 
the kernel is to use the Weinberg rules (20j . which are equivalent to the graph 
techniques in LFD [21 [3]. To calculate the amplitude —Ai, one should put in 
correspondence: 

• to every vertex ~ the factor g, 

• to every intermediate state - the factor 

— -, where Sint = — 

So - Sint + ^0 ^ Xi 

and So is the initial (=final) state energy. 

In our case (three-body bound state): sq = M|, where M3 is mass of the bound 
state three-body system. To every internal line one should put in correspondence 
the factor One should take into account the conservation laws for and 

Xi in any vertex and integrate over all independent variables with the measure 

(27r)3 ■ 

Applying these rules to the graph shown in Fig. [3l we find 

_ f 2rK,2 e{x[)e{x',) d^k\^dx[ 

^ - J s[,,-Mi 2x[ 2x', (27r)3 + i^-^^J + 1-^^^) 



+ E 



2rKi23 e{x[) eix'^) e{x'^) d'^k'i±dx[ d^k'^^dx'^ 

s'l^^-Ml 2x[ 2x'2 2x'3 (M^ (27r)3 



Here x'i + X2 + X3 = 1 in the first term, and x'l -|- X2 + 2:3 = 1 in the last term, 
is the three-body energy in the intermediate state, with two interacting particles 
12, whereas the particle 3 is free (first graph in r.h.s. of Fig. [3]); s'/23 is the three- 
body energy in the intermediate state with all three interacting particles (last 
graph in r.h.s. of Fig. [3]), and F is the LF vertex function. The corresponding 
LF wave function reads 

^(1,2,3) = ^^i^^, 

kl; , +m? kh +m? kn, 
S123 = — + — + — 5 

Xl X2 X3 



where 
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with ki± + k2± + fcsj. = and sum of all x's is 1, as always. -r(l,2,3) means: 
r{l, 2, 3) = r{ki±,xi; k2±,X2; ^3^, X3) and similarly for ■ip{l, 2, 3). 

Now we will use the LF Jacobi momenta k± = ki2±,x = X12 constructed in 
Appendix A. I.e., we make the substitution: 

ki± = k±-xk3±, xi = x(l-X3), 

k2± = -k± - {I - X)k3±, X2 = (1 - x)(l - X3), , . 

A/-j^^ — _L 3^' 1 — V 3/' 

k'2^ = —k'^ — (1 — 2;')fc3_|_, x'2 = (1 — — Xg). 

For (and F) we get -0(1)2,3) = V'(^±) 2;; ^3±) ^s)- Now the variables 
A;_L, x; fc3_L, X3 are independent and they are not constrained by any relations. 
For both x's we have: < x < 1, < X3 < 1. Then S123, Eq. (I5|), becomes 

_ fejx + ^12 fcjx + vr? 

S123 - — ; \ , (7j 

1 - X3 X3 

where 

A;^ + vr? 

S\2 = -7z r 

x(l — X) 

is the effective mass squared of two-body subsystem. The effective three-body 
mass squared S123 has now the form of an effective mass of a two-body system, 
consisting of masses y^si2 and m. 

In terms of Jacobi momenta, the equation for -0 becomes 

(S123 -M|) -0 = / 0ifi23 ^ 



(27r)32x'(l - x') (2^)32x^(1 -x^) 



m? f , d?^' I dx' 



2tt^ J (2^)32x'(l-x')(l-2;3) 
- (231) - (312) . (8) 

Here, instead of K12, given by the Weinberg rules, we have introduced the 
two-body kernel V12 = —-^^Ki2 = y(fc^, x'; fc^, x; M^), which in the non- 
relativistic limit turns into the potential in the Shrodinger equation. The kernel 
y(fc^, x'; fc_L, x; M^) has exactly the two-body form, i.e. the form which enters 
the two-body LF equation. For one-boson exchange it is given in the next section. 



3.1 Scalar LF hound state equation 

In ladder truncation, there is no intrinsic three-body kernel, so we omit the first 
term in Eq. dH) for now. However, we will find the three-body kernel i^i23 due 
to the diagrams of Fig. [5] in Sec. 13.21 below, and calculate its contribution to the 
binding energy perturbatively in the next section. After omitting K123, the LF 
bound state equation, Eq. ([8]), contains three interaction kernels V12, V23, V31. 

For the one-boson exchange model with mass ^ the two-body potential reads 
(see e.g. [2]): 

V{k'^,x';k±,x;M'^) = (9) 

Ana 

f^' + i{l-§fm^ + i{k^-§k'^)' + {x'-x) (^^-M^) 



V.A. Karmanov and P. Maris 



11 



This expression is valid for x < x'. If x > x' , the kernel is obtained from ^ by 
the permutation x <-> x', k± <-> fc^. Here a is the same dimensionless coupling 
constant as introduced in the previous section, a = g'^ /{IGirm?), with g the 
coupling constant in the Hamiltonian H = —gcjp'ip. In the non-relativistic limit 
the kernel, Eq. ([9]), turns into the well-known Yukawa potential 

which in the coordinate space is simply 

V{r) = 



m 
27r' 



r 

Now we introduce the Faddeev components 

V^(l,2,3) = Vi2(l,2,3)+V'23(l,2,3)+V'3i(l,2,3) 

= Vi2(l,2,3)+V'i2(2,3,l)+V'i2(3,l,2), (11) 

such that the component, for instance, ■012(1) 2, 3) satisfies an equation containing 
Vi2 only. This procedure is standard. In this way we obtain the equation for the 
component -012 

(S123 -M|)0i2(fc±,X;fc3X,X3) = (12) 

X 012 {k'i,x';k3±,X3) +ipu (^23±> ^3; M±) a;'i) + 0i2 (4i±, a^si; ^2±) 4) 
where S123 is defined by Eq. ^ and 

^12 = 1 - ^3 M| - ^ . 13 

The three-body mass M| (which we are solving for) enters this equation through 
the factor (S123 — Aff) in l.h.s. and through the value of MI2, Eq. (jl3p . in an 
argument of V . (An equation similar to Eq. (1121) . but for zero range interaction, 
was firstly derived and solved numerically in Refs. [9l 110].) 

The second and third terms in Eq. (fT2]) . namely 0i2 (^23^1 2:23; fc'^, x'j^) and 
■012 (^3i_L) 2:31; ^2^, X2) depend on the LF Jacobi variables (23,1) and (31,2) re- 
spectively, whereas the equation is written in terms of the variables (12,3). The 
variables (23,1) and (31,2) should be expressed though (12,3). These expressions 
are derived in Appendix IA.4[ In particular, the variables A;23_L) ^23' ^i±' -^'i ^'^'^ 
^3i±' 4i' ^2±' ^2 ^^^6 obtained from Eqs. (1A.12|1 and (1A.13P of Appendix I A. 41 by 
the replacement: fc_L — > fc^, x x' . To be specific: 

^, X3k'^ + (1 - x')kss_ ^, _ (l-xO(l-X3) 

^^-^ I - x' + x'x^ ' I- x' + X'X3 ' ^^^^ 

K± = k'^- x'kzL, x'l = x'(l - X3), 
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and 



_ Xsfc^ - X'fcsx , _ X3 



^^-^ + X3 — X'X3 ' x' + X^—x'x^^ ^-j^g-j 

= -K - (1 - 4 = (1 - x')(l - xa). 

With this notation (and omitting prime), the full wave function obtains the form: 

^p{k±,x; k3^,X3) = ^pi2{k±,x; fcg^jXs) (16) 

+ Tpl2{k23±,X23; ki^,Xi) + ljji2{k3i±,X3i; k2±,X2) ■ 

It is normalized as 

, 2 d'^k±dx (fk^^dxj, 

J iV'lfc^'^i^a^'^s)! (2.)32x(l-x) (2^)32x3(1 -X3) = ' " ^'^^ 



3.2 Contribution from three-body kernel 

Naively, one might not expect any three-body kernels in the ladder truncation. 
However, in LFD, there are three-body kernels, as already mentioned in the in- 
troduction, due to the (LF) time-ordering of the diagrams. To find the correction 
to the binding energy from the kernel K123, we represent equation for the full 
wave function symbolically in the form 

(S123-M|)V' = iki2 + k23 + K3l)i^ + Ki23i^ . 

Substituting M| ^ M| + Z\M|, we find 

AMi = -i^kusi'), 



or, explicitly 



^ J (27r)3 (2^)3 (27r)3 (27r)3 ^ ^ 



X 



^ dx3 dx dx'n dx' 



2x3(1 -X3) 2x(l-x) 2x^(1 -x^) 2x'(l-x') 
x^{k^,x; k3±,X3) Ki23{k_L,x; k3±,X3; fc^, x'; ^3^, Xg) ■0(^1, x'; ^3^, X3) . 

The integrand depends on four 2-dimensional vectors k±, k3±, k'^, feg^ and, 
hence, on three relative angles between them. Therefore Eq. (llSp is an 11- 
dimensional integral (the last, fourth angle integration, is reduced to a factor 
2tt). 

The kernel -ftri23 is determined by 54 graphs: 48 of them contain two exchange- 
bosons in flight, whereas another 6 correspond to the intermediate production of a 
pair of constituent particle-antiparticles. It is sufficient to calculate only the nine 
graphs shown in Fig. [2j The other contributions are obtained by permutations 
and can be taken into account by multiplying the results of these nine graphs by 
a factor 3!=6. 
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Here we show explicitly how to calculate the contribution of the first graph 
Fig. [2] only. All other contributions are calculated in a similar fashion. Following 
the Weinberg rules [20], we get 



K 



(1) 

123 



287r2m4a2 



0{xi 



{x'l + x'2 - Xi) 9{X3 



X'n 



(^1 



Mi) 



(^2 



Mi) 



[S3 



^i) {xi-x[) {x[+X2-xi) {x3-x'^) 



(19) 

where Si are the energies in intermediate states shown in Fig. [2l The line 1, 
for instance, is associated with the momenta fci_L,xi and contributes in si the 

term -^-^^ — , and analogously for other lines. We sum over all the lines in given 



Xl 



intermediate state. In this way we find for si: 



si 



Xl 



X2 



X3 



+ 



k'lj_ + 



Other intermediate state energies can be found similarly. Transforming this ex- 
pression to the LF Jacobi variables by the substitution we find: 



si 



3± 



+ 



fc^ +x(l 



x)k |_|_ 



m 



X3 



x{l - x)(l - X3) 



+ 



The final expression for the kernel ^^123 obtained by substitution of this 
expression for si (and analogously, for ^2,53) in Eq. (fT9|) and replacing also the 
variables x's in Eq. (jl9p by the Jacobi ones according to Eqs. ([6]). The 0-functions 
impose the following constrains on the Jacobi's x's: 



< Xs < 1, < X < 1, < Xo < X3, 



and 



< x' < 



x{l - X3) 
1-x', 



Other contributions are found similarly. All of them have the form (jl9p with 
the replacement of si, 82,83 by the corresponding twelve intermediate energies 
•Si — 812, as shown in the various graphs in Fig. [2l and also x's in the 0-functions 
and in denominator. Taking the sum of all the contributions and substituting it 
in Eq. ([18]) . we find the perturbative correction to the binding energy resulting 
from the three-body forces: 



AMi 



+ 



-6 



(fk3± (Pk± (fk'^^ (fk'^ 
~{2^ (2^ "(2vrF 
^3 dx'o 



2x3(1 



V^(KS + i^i|)V''dx' 



2x'(l - x') 



+ 



24(1 - 



dX3 



dx 



2x3(1 -X3) 70 2x(l 



v.(i^S + ^iS 



dx' 



l-x' 



2x'(l 
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2x'{l - x') h(i-^a) 2x'{l - x') 

l-x'^ J 

dx's tpK^^^jj' dx' 

+ io 2x'^{l-x'^) Jl^ 2x'{l-x') 

l-x^ 

The three-body wave function ip is expressed through the Faddeev component 
■012 by Eq. p^ . whereas ■012 is found from Eq. (fT2|) : ^p and ^p' depend on the 
variables k±,x; k^±,x^ and k'j_,x'; k'^j_,x'^ respectively. 

4 Numerical results 

The LF wave functions in Eq. (|12p depend on five independent scalar variable^. 
The wavefunction tpi2ik±, x; k^x, xs) on the l.h.s. of Eq. (fT2]) depends on the 
external variables k±, k3±, x, x^, and z = cos /3 = k± • k3_i/{k_ik3±). The argu- 
ments of the wave functions ^jJ''s in the integrand on the r.h.s. of Eq. (jl2p can be 
expressed through a combination of these five external variables and the three 
integration variables k'j_, x', and cos = k± ■ k'j_/{k±k'j_). Thus the structure of 
the LF bound state equation, Eq. (|12p . is very similar to that of the BS equation, 
Eq. and we can solve it in a similar fashion. That is, we use a transformation 
to map the radial variable onto a finite interval, discretize all five variables, and 
perform all integrations numerically (three in the case of LED, four in the case 
of the BS equation). 

4.1 Difference between BS approach, and LFD 

For small binding energies, our results for the three-body binding energies = 
3m — M3 are shown in Fig. [5] (the units are set by m = 1). For comparison, we 
also include the non-relativistic result obtained by solving the non-relativistic 
Schodinger equation (solved via corresponding non-relativistic Faddeev equation 
in the momentum space with the kernel Eq. (jlOp ). Clearly, all three calculations 
give different results. 

However, the relativistic binding energies, calculated using the BS equation 
and in LFD, are approaching to each other as the binding energy goes to zero. 
Unfortunately, the numerics become less accurate in that limit, because both 
the interaction kernels and the wave functions we are solving for become more 
peaked. Nevertheless, for /x = 0.5, within numerical accuracy, they appear to 
have the same critical coupling, defined as the coupling constant a at which 
the binding energy becomes zero. In contrast, the non-relativistic binding energy 
behaves quite differently, and vanishes at a much smaller value of a. 

This behavior is similar to that which was found in the two-body system 
In the latter case, the decreasingly small binding energies, calculated via 
relativistic (both BS equation and LFD) and non-relativistic equations, become 
close to each other only in the limit /i — > 0. 

^We remind that three-body BS amplitude also depends on five scalar variables. This coinci- 
dence takes place for two- and three-body systems only. For n-body system with n > 4 LF 
wave function and BS amplitude depend on different numbers of scalar variables. 
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dimensionless coupling constant a 



Figure 5. (Color online) Three-body binding energy in the limit when it tends to zero, for 
jj, = 0.01 and fi = 0.5 (units set by the constituent mass: m = 1). Solid curves are the BS 
result; dotted curves are the LFD results; dashed curves are the solution of the non-relativistic 
Schrodinger equation. 



Indeed, for the value /i = 0.01 in Fig. [5l the relativistic and non-relativistic 
results approach each other as all the binding energies decrease down to Ef, ~ 
0.01, and they seem to extrapolate to a ~ as the point at which the binding 
energy becomes zero. We did not calculate for smaller Eh and /i because of 
increasing numerical errors. 

More intriguing is the rather large difference we see in this figure between the 
results from the BS equation and those from LFD for nonzero binding energy Eb. 
The LF equation in ladder truncation significantly underbinds the three-body 
system relative to the ladder BS equation. Though at small binding, both BS 
and LFD binding energies tend to each other and to zero, the relative difference 
between them remains significant, and is roughly a factor of two. We remind the 
reader that for two-body system these binding energies approximately coincide 
with each other [3]. 

A closer look at our results reveals that the difference between the three-body 
masses squared in the BS approach and in LFD, j^p—M^ ^g, is approximately 
proportional to itself. In Fig[6l we show the difference 

AM^ = mIlp - Miss 
as function of 9m^ — M|^^. It is a remarkably linear function over the entire 
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accessible interval of (remember: 9m^ — M| = means zero binding 

energy, and M| = is the maximal binding energy); furthermore, it depends 
only weakly on /i, the mass of the exchange boson. 



4.2 Differences largely explained by 3-body force 



The next question is: Can this difference between the BS approach and LFD be 
explained by the diagrams of Fig. [2]? In order to answer this question, we per- 
formed a perturbative calculation of the correction to due to these three- 
body forces, using Eq. ()20|) . In Fig. [7| we show our results for the contributions 
both from the diagrams with two exchange-bosons in flight (first eight graphs 
in Fig. [2]) and from the diagrams with creation and annihilation of an interme- 
diate particle-antiparticle pair (ninth graph in in Fig. [2]) . Since this calculation 
involves the numerical evaluation of an 11-dimensional integral, see Eq. (Iisp . the 
results are not highly accurate; nevertheless, they are quite illustrative. 

As we can see in the top panel of Fig. [71 the contribution of two exchange- 
bosons in flight to AM^ is an almost linear function of Mg. Furthermore, it 
depends only very weakly on the value of // in the wide interval 0.01 < /i < 3.0. 

The contribution of the pair creation term only (ninth graph in Fig. [2|) is 
shown in the bottom panel of Fig. [71 Again, it is almost linear in M|, but, in 




9 m - Ml 



Figure 6. (Color online) The difference M^j^p — vs. the deviation of M^ gg from the 

free value: 9m^ — M^gg, for different values of the exchange mass ^. 




1 2 3 4 5 

3, LF 

Figure 7. (Color online) Correction ZiMf, Eq. ([201, vs. the difference 9m^ - M^^^pj^ for 
different values of the exchange mass ^. Top: Contributions from the first eight graphs in Fig. [2] 
(without pair-creation); Bottom: Contributions from pair-creation only, the last diagram in 
Fig. [2] Notice the difference in the vertical scale. 
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contrast to two exchange-bosons in flight, it strongly depends on ji. At ^ = 3.0 
this contribution is about 20% relative to two exchange-bosons in flight but it 
becomes negligible for /x < 0.15. 

The behavior of these both contributions is consistent with the behavior of 
the difference of the masses squared in the BS and LFD approaches, Fig. [H This 
linearly increasing (in absolute value) negative correction from effective three- 
body forces in LFD eliminates most, if not all, of the linearly increasing positive 
difference between BS and LFD results. 

The three-body bound state mass squared M| for the exchange masses ^ = 
0.01, 0.5 and 1.5 are shown in Fig. [HI Like in Fig. [5l the solid and dotted curves 
are our results from the three-body BS equation in ladder truncation, Eq. (j4]), 
and our results from the LF 3-body bound state equation with the two-body one- 
boson exchange kernel, Eq. Q, respectively. As in Fig.O the BS and LFD results 
differ significantly from each other and the difference increases with increase of 
the binding energy (and, hence, with the coupling constant). 

The crosses on the dotted curves indicate the points where we evaluated 
the perturbative correction due the 3-body forces depicted in Fig. [2l The dia- 
monds with error bars in Fig. [8] indicate the sum M| + AM^ , where AM^ 
is the correction due from the diagrams of Fig. [2l calculated perturbatively by 
Eq. ()20p , and displayed in Fig. [71 The errors indicate the numerical uncertainty, 
mainly due to the numerical evaluation of an 11-dimensional integral. This cor- 
rection does indeed shift the LFD three-body mass squared so that it comes in 
reasonably good agreement with the BS results. This is the effect of relativis- 
tic three-body forces, dominated here by two exchange-bosons in flight, with a 
small contribution from the creation (and annihilation) of a pair of constituent 
particle- antiparticles. 

One should of course keep in mind that here we only present a perturbative 
estimate of the correction due to these effective 3-body forces in LFD. Note 
that the perturbative calculation is unexpectedly good even for large correction 
compared to the value of itself, especially for small values fi of the mass of the 
exchange particle. The agreement becomes less pronounced as the exchange mass 
increases, once the exchange mass becomes larger than the constituent mass. 
E.g. if the exchange mass is three times the constituent mass, the perturbative 
correction explains about 70% of the difference in M| between the BS approach 
and LFD. 

Since the corrections are substantial, one could raise the question: How is a 
perturbative calculation reliable in this context? To answer this question, one 
should solve equation ([8]) beyond the perturbative framework, which is numeri- 
cally much more demanding and requires signiflcantly more work. Nevertheless, 
we have ample evidence that the contributions of effective 3-body forces in LFD 
are substantial in the three-body systems, and can explain most of the difference 
between the results found in the BS approach and in LFD. 



V.A. Karmanov and P. Maris 19 




dimensionless coupling constant a 




dimensionless coupling constant a 



Figure 8. (Color online) Three-body bound state mass squared Af| vs. coupling constant a 
for exchange masses /i — 0.01 (top), 0.5 and 1.5 (bottom). The units are set by the constituent 
mass: m = 1. 
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5 Conclusion 

We have solved, for the first time, the three-body BS and LFD bound state 
equations with a one-boson exchange kernel and found a difference in the cor- 
responding binding energies. This difference is absent in the two-body BS and 
LFD equations [4J. 

When the binding energies tend to zero, the absolute difference between BS 
and LFD results decreases too. However, the difference of the binding energies 
relative to the binding energy itself remains roughly constant: LFD in ladder 
truncation underbinds approximately by a factor of two, compared to the lad- 
der BS equation. The relative difference depends only weakly on the exchange- 
boson mass /i. At the same time, the difference between the relativistic BS (or 
LFD) binding energies on one hand, and the binding energy found from a purely 
non-relativistic approach on the other hand, is much larger than the difference 
between the BS approach and LFD. In particular the critical coupling for the 
onset of bound states in a three-body system appears to be the same in the 
BS approach and in LFD, but in a non-relativistic calculation this critical cou- 
pling is significantly smaller. This difference decreases with the decrease of the 
exchange-boson mass; in the limit of a massless exchange boson, this critical 
coupling tends to zero for all three approaches. The latter is quite similar to 
two-body calculations, where one also finds a significant difference between a 
purely non-relativistic calculation on the one hand, and the BS equation or LFD 
on the other hand. 

Most of the difference between BS and LFD three-body binding energies can 
be attributed to effective 3-body forces. Despite the fact that both in the BS 
approach and in LFD we use an one-boson exchange kernel, the actual kernels 
in the three-body BS and LFD bound state equations differ from each other by 
contributions of irreducible LFD graphs with two exchange-bosons in flight and 
by pair creation. Both types of graphs are implicitly included in the ladder kernel 
of the BS equation (this becomes obvious if its iterations are transformed in a 
set of the time-ordered diagrams), but they are not included in the one-boson 
exchange kernel in LFD. In the latter case they should be taken into account 
separately, as an extra contribution from a 3-body force. After taking them into 
account perturbatively, the LFD binding energy becomes approximately equal 
to the BS results, at least for exchange masses smaller than the constituent mass 
(fi < m). This is a clear manifestation of three-body forces, which in our model 
are unambiguously given by the set of graphs Fig. [2] and explicitly calculated. 

Three-body forces of this origin contribute only to the kernel of 3-dimensional 
bound state equations, like the LFD one. They can also be taken into account as a 
relativistic correction in the framework of the Schrodinger three-body equation 
(together with other relativistic corrections). They are true three-body forces 
in a 3-dimensional framework. However, in the framework of the BS equation, 
using 4-dimensional integral equations and explicitly covariant variables without 
any 3-dimensional reduction, these effects are automatically incorporated via 
iterations of the two-body ladder kernel. They do not manifest themselves as 
three-body forces in the BS approach. Therefore, the interpretation of forces as 
two- or three-body forces, depends on the framework with which these forces are 
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associated. 

Note that three-body forces, resulting from excitations of intermediate isobar 
states, do not depend on the type of equation: BS, LFD or Schrodinger equation. 
They remain three-body forces in any of these equations. However, they may turn 
into two-body forces if one considers the intermediate isobars in the framework 
of a system of coupled channel equations, including transitions N <^ A. Both 
types of these three-body forces could be taken into account as a correction to a 
non-relativistic treatment of nuclei. 

In addition, there can be intrinsic three-body forces, such as the six-nucleon 
contact term NNN — > NNN, generated by chiral perturbation theory for the 
effective nucleon interaction at the N2L0 level [21j . or, in QCD, three-body 
forces due to the triple gluon vertex. However, these intrinsic three-body forces 
are beyond the skope of the present paper. 

We would like to emphasize that although we find a correction to the LFD 
kernel (three-body forces) which explains most of the difference between the 
LFD binding energy and the BS binding energy, this does not mean that LFD 
is considered as an approximate method to solve the BS equation. Both the BS 
approach and LFD have their advantages and disadvantages. Without trunca- 
tions, they should give the same results for physical quantities such as the bound 
state mass. However, truncations can lead to differences between LFD and the 
BS approach. A detailed comparison between results obtained in LFD and in the 
BS approach, such as the work presented here, provides deeper insight in the lim- 
itations and deficiencies of the truncation. Here, it revealed that for three-body 
bound states in the ladder truncation in LFD there are significant contributions 
to the binding energy coming from effective three-body forces related to higher 
Fock sectors. 

Finally, the calculations presented here have been done for spinless systems: 
a scalar bound state, composed of three scalar constituents, interacting via ex- 
change of a scalar boson. Physical systems such as the triton and ^He nuclei, 
or baryons as bound states of three quarks, consist of fermions, interacting via 
(pseudo-)scalar and/or vector bosons. We expect that at least qualitatively our 
conclusions hold for those systems as well: For three-body systems, the ladder 
truncation of the Bethe-Salpeter equation and the ladder truncation in light- 
front dynamics are not equivalent, and are likely to give significantly different 
results, due to the absence of diagrams like those in Fig. [2] in the light-front 
ladder truncation. 
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Appendix A: LF Jacobi variables 

The three-body LF wave function depends on the variables fcix, fc2±, fcsx and xi, X2,xz- Because 
of conservation of momenta, they are not independent, but satisfy the relations: feix + fe2x + 

fcsx = Q, Xi + X2 + X-i = 1. 

For any sub-system (a/3) we define corresponding LF Jacobi variables, which are indepen- 
dent from each other. 



A.l Subsystem 12 

For subsystem 12 we introduce the Jacobi momenta as follows: 

X1+X2 ' (A 1) 

fei2± = fcix - a;i2(feix + fe2x) • 

The construction of xi2 agrees with general definition of variable x: for subsystem 12 it is the 
ratio Xi2 — j^^, where pi2+ = Pi+ + P2+ is the plus-component of the four-momentum of 
the subsystem 12. Dividing numerator and denominator by the plus-component of the total 
three-body momentum p+ = po + Pz — pi+ + P2+ + P3+ , we find for xi2 expression HA.1|I in 
terms of the variables xi,X2 defined for three-body system. 
One can express them as follows 

^ 1^ ' (A 2) 

fel2± = fcl± + T^fe3± ^ ' ' 

and vice versa: 

11 = a:;i2(l - 2:3), feix = fei2± - a::i2fc3x- (^-3) 

The Faddeev component -012 becomes a function of 

V'12 = ^/'i2(fci2±, 2:12; feax, 3:3). 

In contrast to the variables xi,a;2,a;3, constrained by xi + X2 + X3 — 1, both variables xi2,X3 
vary independently in the interval from to 1. 

The three-body integration measure, with integration over two-body subsystem, is trans- 
formed as: 

r(2)/,/ , ,/ , , \xr ' < '1 .d^k'ij^dx'id^k'2±dx'2 

0^ '(feix + fe2± + fe3±)0(a;i + X2+X3 - 1) — 



(5'^'(fe'l2± + k2i^)5{x[2 + X21 - 1)- 



23^1 x^ 

1 d k^2j^dxi2d k2\A_dx2\ 



(1 X-i) 2x-^22;2l 

d^k'i2^dx'i2 



(I-X3) 2x'^^(l - x'^^) 

In the non-relativistic limit X12 ~ 1/2, xz ~ 1/3, and we obtain usual Jacobi coordinates: 



fel2X ~ feix - ^(fcl± + fe2x) = ^(felX - fe2±), 



1 2 

fe3X ~ feix - -(fclX + fe2X + fesx) = - 



'^i-L - ^C^ax + fesx) 



(A.4) 



A.2 Subsystem 23 



Corresponding formulas are obtained from the relations for the 12-subsystem by the cyclic 
permutation 123 231. Then the Jacobi variables are defined as follows: 

^23 = , 

fe23X = fe2x - a;23(fc2X + fesx) • 
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One can express them as functions of the three-body variables 
and vice versa: 



a:2 = 2:23(1 - a::i), fc2± = fe23± - a;23fei±- (A.7) 
The Faddeev component 'i/'23 becomes a function of 

^"23 = V'23 (fc23± ,X23;ki±,Xl). 

A. 3 Sub-system 31 

Corresponding formulas are obtained from the relations for the 23-subsystem by the cyclic 
permutation 123 231 or from the relations for the 12-subsystem by the cyclic permutation 
123 — > 312. Then the Jacobi variables are defined as follows: 



X3+X1 ' 

fe3i± = fcsx - a;3i(fc3i + fcix) • 

One can also introduce 2:13 = 1 — X31, fci3x ~ — fcsix- 

One can express them as functions of the three-body variables 

k3i± = fc3x -I- Y^^k2± 

and vice versa: 



(A.8) 



(A.9) 



2:3 = 2:31(1 - K2), fc3X = fe3ix - 2:3ifc2X- (A. 10) 

The Faddeev component tpsi becomes a function of 

'031 = ^^31 (fc31± , 2:31 ; fe2X , 2:2) . 

A. 4 Relation between different Jacobi coordinates and cyclic permutations 

We express the coordinates for the subsystems 23 and 31 in terms of 12: 2:12, fei2±, 2:3, fe3x. For 
that we substitute in the above formulas the Eqs. (|A.3|I and use the relations: 



(A.U) 



Then we find: 



2:1 = 2:12(1 - 2:3), feix = fei2X - 2:i2fc3x, 

2:2 = 1 - 2:12(1 — 2:3) — 2:3, fe2X = -(fel2± - 2:i2fe3x) — fe3i 



(1 - 2:i2)(l - 2:3) 2:3fei2X + (1 - 2:i2)fc3 

2:23 = — : , fc23± ~ 



and 



1 — 2:12 + 2:122:3 1 — 2:12 -I- 2:122:3 

2:1=2:12(1-2:3), feix = fel2± - 2:i2fe3±- 



2:3 , 2:3fci2X - 2:i2fe3 

2:31 = ; , K31X = 



(A.12) 



2:12 + 2:3 — 2:122:3 2:12-1-2:3-2:122:3 

(1 - 2:i2)(l — 2:3), A;2X = -fel2± — (1 — 2:i2)fe3J 



(A.13) 



After cyclic permutation 123 231 the variables fei2x,2:i2 turn into fc23x,2:23 defined by 
Eqs. I|A.12|I , and similarly for other variables and permutations. In this way one obtains Eq. (|16p 
for total wave function. Note that the transformations of the 2:-variables are non-linear, since 
they are defined through 2:1,2:2,2:3 by the non- linear formulas (|A.ip . (|A.5p . (|A.8p . 
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